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The spin dynamics of electron doped Mott insulators on a triangular lattice is studied based on 
the t-J model. It is found that the particularly universal behaviors of integrated dynamical spin 
structure factor seen in the doped Mott insulators on a square lattice, are absent in the doped 
Mott insulators on a triangular lattice, indicating the presence of the normal state gap. As a result, 
the spin-lattice relaxation rate l/Ti divided by T reduces with decreasing temperatures in the 
temperature region above 0.2Jf5i50K, then follows a Curie-Weiss-like behavior at the temperature 
less than 50K, in qualitative agreement with experimental observations. 



PACS numbers: 71.70.Ej, 73.20.At, 74.20.-z 



INTRODUCTION 



In the past few years, the discovery of superconduc- 
tivity in layered cobalt oxides Naa;Co02 ■ yH20^ has at- 
tracted considerable attention among the condensed mat- 
ter physics society. Although the superconducting tran- 
sition temperature (Tc w 4.5K)^ is very low, the layered 
electronic structure and the dome-shaped behavior of Tc 
suggest this hydrated compound may be another doped 
Mott insulator—!^. In contrast to the cuprate supercon- 
ductors with a square lattice of Cu02 planes, the par- 
ent compound Naa;Co02 is composed of two-dimensional 
(2D) C0O2 layer with triangular geometry. This ge- 
ometry frustration results in some novel electronic and 
magnetic phases, for instance, the Anderson's resonat- 
ing valence state"' and the strong topological frustration 
phases^'^. With the doping variation, Naa;Co02 displays 
rich and complicated phase diagram''. In particular, su- 
perconductivity appears around x = 0.3 when water 
molecules are intercalated between the C0O2 layers^i^. 

Although extensive researches on this compound have 
been done in both experimental and theoretical sides, 
the doping dependent magnetic properties still remain 
unclear. The local spin density approach^ predicted fer- 
romagnetic coupling within 2D C0O2 layer for nearly all 
electron doping 0.3 < x < 0.7. In fact, up to now, only 
at large doping neutron scattering studies and time-of- 
flight experiments on Nao.75Co02 and Nao.82Co02^iiS 
respectively, present low-energy fluctuations characteris- 
tic of A-type antiferromagnetism: antiferromagnetically 
coupled ferromagnetic layers. However, with the decreas- 
ing of Na content, a series of nuclear magnetic resonance 
(NMR)^^ and nuclear quadrupole resonance (NQR)-^ as 
well as Knight shift measurementfiS show the presence of 
strong 2D antiferromagnetic (AF) spin correlations and 
suggest the possibility of the spin-singlet superconduc- 
tivity. That is to say the 2D AF spin correlations may 
be realized in the region of small x value for Na2,Co02. 
Starting from ab-inito band structure calculations, Kor- 
shunov and co-workers have estimated a certain critical 
concentration Xm, below which the magnetic suscepti- 
bility within the C0O2 plane shows a tendency towards 



AF fluctuations^. The recent correlation effects on the 
doped Naa;Co02 via a cellular cluster approach has recon- 
firmed and explained such in-plane magnetic transition 
from FM tendencies towards to AF— . 

Recently, the surface-sensitive angle resolved photo 
emission spectroscopy (ARPES) experiments^^"— ob- 
serve a doping-dependent evolution of the fermi surface, 
which is found to be centered around the F point with 
a mostly aig characterister and no sign of six pockets 
resulting from the band predicted by local density ap- 
proximation calculation^ for a wide of Na concentrations. 
This indicates the importance of the electronic correc- 
tions in the doped Naa:Co02. Therefore, as discussed 
by many researcher o^i^'^^i^"' , an effective single band low- 
energy model suffices to capture the essential physics of 
electron-doped cobaltate Naa;Co02. With this considera- 
tion in mind, in this paper we also take a single band t-J 
model as our starting point to discuss its unusual dynam- 
ical spin response. Since the AF fluctuations within the 
C0O2 planes may be realized in the low electron doping 
range, we are focusing on such region in our calculations 
below. We find that the particularly universal behaviors 
of integrated dynamical spin structure factor seen in the 
doped cuprates, are absent in the doped cobaltates, indi- 
cating the presence of the normal state gap. We also find 
that the spin-lattice relaxation rate 1/Ti divided by T re- 
duces with decreasing temperatures in the temperature 
region above 0.2Jw50K, then follows a Curie- Weiss-like 
behavior at the temperature less than 50K originating 
from the strong AF fiuctuations, in qualitative agree- 
ment with experimental observation !] ^ ^ 1 . Our results 
also show that AF fiuctuations within the C0O2 plane 
gradually weaken with increasing electron doping. 



II. THE t-J MODEL AND FERMION-SPIN 
THEORY 



It has been argued that the essential physics of the 
doped C0O2 plane is contained in the t-J model on a 
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triangular lattice^, 

irj(7 ia 
+ J Sj ■ Sj-i-jj, (1) 
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where the summation is over all sites i, and for each i, 
over its nearest- neighbor fj, C}^ (Cia) is the electron cre- 
ation (annihilation) operator, Si = CjaCi/2 is the spin 
operator with a = {<Jx,<^y,<^z) as the Pauli matrices, /i 
is the chemical potential, and the projection operator P 
removes zero occupancy, i.e., '^^Cj^Cia > 1. In the 
past fifteen years, some useful methods have been pro- 
posed to treat the no double occupancy local constraint 
in hole doped cuprates. In particular, a fermion-spin 
theory based on the partial charge-spin separation has 
been developed to study the physical properties of doped 
cuprate a^^i^^ , where the no double occupancy local con- 
straint can be treated properly in analytical calculations. 
To apply this theory in the electron doped cobaltates, the 
t-J model (1) can be rewritten in terms of a particle- hole 
transformation C,, 



H 



Si- a as. 



supplemented by the local constraint fiafi<^ — 1 
to remove double occupancy, where for convenience, we 
have set t = —t^ < 0, fj^ (fia) is the hole creation 
(annihilation) operator. Then the hole operators can 

be decoupled as, fif = al^S~ and /ij. = al^S^ , in the 
charge-spin separation fermion-spin theory22i2^, where 
the spinful fermion operator aia = e~'*'''ai describes 
the charge degree of freedom together with some ef- 
fects of the spin configuration rearrangements due to 
the presence of the doped electron itself (charge car- 
rier), while the spin operator Si describes the spin de- 
gree of freedom (spin), then the single occupancy local 
constraint, 'ZafLf^'^ = S+a^^al^S;- + 5^0,4,0 = 

aial{S^ S~ + S^) = 1 — ajoi < 1, is satisfied in analyt- 
ical calculations. In this charge-spin separation fermion- 
spin representation, the low-energy behavior of the t-J 
model (2) can be expressed as^, 



^" ■^ is the 
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with Jeff = (1 — x)^ J, and x 
electron doping concentration. In this case, the mag- 
netic energy (J) term in the t-J model is only to form 
an adequate spin configuration, while the kinetic energy 
{t) term has been transferred to the charge carrier-spin 
interaction, which dominates the essential physics. 



In the framework of the charge-spin separation, the ba- 
sic low-energy excitations are charge carriers and spins. 
It has been shown that the charge dynamics is mainly 
governed by the scattering from the charge carriers due 
to the spin fluctuation^^, while the spin fluctuations cou- 
ple only to the spin and therefore no composition law is 
required in discussing the spin dynamics^Si, but the effect 
of charge carriers is still considered through the charge 
carrier's order parameter entering the spin propaga- 
tor. In this case, the spin dynamics of the doped square 
antiferromagnet has been discussed^^ by considering the 
spin fluctuation around the mean-field solution, where 
the spin part is treated by the loop expansion to the sec- 
ond order. Following their discussions, we can obtain the 
dynamical spin structure for the electron doped cobal- 
tates as 

/•OO 

S{k,uj) = Re / dte''"^'~*">D{k,t-t') 
Jo 

= -2[l + nB{uj)]lmD{k,uj), (4) 

respectively, where the full spin Green's function, 
D-^{k,uj) = D(°)-i(fc,a;) - T,s{k,uj), with the mean- 
field (MF) spin Green's functions D(°)-i(fc,w) = (w^ _ 
LoD/Bk, and the second-order spin self-energies from the 
charge carrier pair bubble, 
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S.p' ~ ^k+p ^ + £.p+p' ^ Cp' + '-^k+p 



where Z is the number of the nearest neigh- 
bor sites, Bk = X[2x^{€-fk - 1) + xilk - e)], 
A = 2ZJoff, e = 1 + 2*0/ Jeff', 7fc = [cosk^ -\- 
2cos{kx/2) cos{\/3ky/2)]/3, the spin correlation tunc- 
tion X" = (^fSf+f,>, Fl'\k,p,p') = nF(W)[l - 
nF{^p')]-nBiuJk+p)[nF{Cp')-nFi^p-\-p')], F^''\k,p,p') = 

"f(Cp+p')[1 - "-f(Cp')] + [1 + "•B(wfe+p)][nF(Cp') - 

nF(Cp+p')]i nsi^jp) and nF(Cp) are the boson and fermion 
distribution functions, respectively. The MF charge car- 
rier and spin spectra are given by ^fe = Ztxjk ~ M; and 
Lol = Ai(7fc)2 -I- A2jk + A3, respectively, with Ai = 
aeA2(ex^-hx/2), A2 = -eA2[a(x^ + ex/2) + (aC^ + (1 - 
a)/(4Z) - aex/(2^)) + (aC -h (1 - a)/{2Z) - ax' 12)1 2^ 
Aj, = \\oiC' + (1 - a)/(4Z) - aex/(2^) + <?{(xC + 
(1 — Q;)/(2Z) — ax^/2)/2], and the spin correlation func- 
tion X - (^+5r+^), C = (l/^2)E,W-(^if,^,r+^-,>, and 
C' = (l/^2)E^,,'(^^+f,^r+,')- o'^dcr not to violate 
the sum rule of the correlation function (S^S^) = 1/2 
in the case without AF long-range order, the important 
decoupling parameter a has been introduced in the MF 
calculation2Si21^ which can be regarded as a vertex cor- 
rection. All the above MF order parameters, decoupling 
parameter a, and chemical potential /i are determined by 
the self-consistent calculation^i. 
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III. NUMERICAL RESULTS AND 
DISCUSSIONS 
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FIG. 1: The integrated dynamical spin structure factor at the 
doping (a) x = 0.33 and (b) x = 0.30 with the temperature 
T = 0.3 J (dashed-dotted hne), T = 0.4 J (dotted hne), and 
T = 0.6J (solid line) for the parameter t/J — —2.5. Inset: 
the numerical simulation taken from Ref. [27]. 



120 



100 - 




0.4 0.3 0.2 0.1 0.0 0.1 0.2 0.3 0.4 

a)/|t| 

FIG. 2: The integrated dynamical spin structure factor at the 
doping (a) x = 0.33 and (b) x = 0.30 with the temperature 
T = 0.7 J (dashed-dotted line), T = 0.8 J (dotted line), and 
T = 0.9J (solid line) for the parameter t/J = —2.5. 

For understanding the dynamic spin response of the 
electron doped cobaltates, we have calculated the inte- 
grated dynamical spin structure factor, which can be ex- 
pressed as 



The numerical results of the integrated dynamical spin 
structure factor at doping (a) x — 0.33 and (b) x = 0.30 
with the temperature T = 0.3 J (dashed-dotted line), 
T = 0.4J (dotted line), and T = 0.6J (solid hne) for the 
parameter t/3— -2.5 are plotted in Fig. 1. It is shown 
that the integrated dynamical spin structure factor de- 
creases with increasing energies for ui < 0.3t, and almost 
constant for ui > 0.3t. which is similar to the case of 
doped cupratea2^>^ and reflects the common features of 
the doped Mott insulators. However, the temperature de- 
pendence of the integrated dynamical spin structure fac- 
tor, which is essentially temperature independent at high 
energy in a wide temperature regime, presents a strong 
deviation at the low-energy range, so that the shape of 
the integrated dynamical spin structure factor does not 
follow a universal behavior, in contrast to the case of the 
doped cuprates (seen the insert in Fig. Ib)^^. It has been 
discussed that the universal behavior of the integrated 
dynamical spin structure factor in the doped cuprates is 
due to the absence of the normal state gap in the electron 
density of states^^. However, in the doped cobaltates, the 
normal state gap open a^^i^^'^" , which is mainly induced 
by the frustrated spin. Since in the charge-spin separa- 
tion framework the spin dynamics is dominated by the 
scattering of spins, which are strongly renormalized be- 
cause of the strong interactions with fluctuations of sur- 
rounding charge carrier excitations. The frustrated spin 
moves in the background of the charge carriers, and the 
cloud of distorted charge carriers background is to follow 
the frustrated spins, which leads to the anomalous spin 
dynamics in the doped cobaltates. In other words, the 
origin of the absence of the particularly universal behav- 
ior of the integrated dynamical spin structure factor in 
the doped cobaltates at low energy originates from the 
normal-state gap due to the magnetic frustration. In or- 
der to show this normal state gap clearly, we calculate 
the S{uj,T) at higher temperature as shown in Fig. 2, 
where the parameters are the same with the Fig. 1 ex- 
cept for the temperature T = 0.7 J (dashed-dotted line), 
T = 0.8J (dotted line), and T = 0.9J (solid line). In Fig. 
2, the temperature dependent integrated dynamical spin 
structure factor nearly degenerates and doesn't show a 
deviation at low energy seen in Fig. 1. This universal be- 
havior of the integrated dynamical spin structure factor 
exists for the whole temperature range is similar to the 
case of doped cuprates (seen the insert in Fig. Ib) ^^'^^ . 
From Fig. 1 and Fig. 2, we estimate the normal state 
gap opens at temperature Tgap — 0.7 J. According to our 
estimation J ~ 250K—, we get the Tgap — 160K, very 
close to the experimentsi^ . 

Now we turn our attention to discuss the spin-lattice 
relaxation rate l/Ti. Theoretically, the spin contribu- 
tions to the spin-lattice relaxation rate 1/Ti divided by 
T may be written using the imaginary part of the Co 
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dynamical electron spin-susceptibility x"(k, w) as 



1 



lim 



\A{k)\ 



2X"(k,^) 



(7) 



where x"{k,io) = [1 + risii^)]^^ S{k,uj), g is the g factor, 
Hb is the Bohr magneton, and A{k) — J^i Aie^'^^\ where 
Ai is the hyperfine coupling between the electron spin and 
the nuclear spin having gyromagnetic ratio ji^^ and is to 
be approximated to a constant in our calculation. The 
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FIG. 3: (color online) The spin-lattice relaxation rate l/Ti 
divided by T of the electron doped cobaltates as a function of 
temperature at doping x = 0.30 (solid line), x — 0.33 (dashed- 
dotted line), and x = 0.40 (dotted line) for parameter t/J — 
—2.5, where we choose 27^ fc_Bj4i/(;^/i^=constant. Inset: the 
experimental result taken from Ref. [10]. 



spin-lattice relaxation rate I/Ti divided by T in Eq.(7) 
has been evaluated numerically and the results for t/J= 
-2.5 at the doping x=0.30 (solid line), x=0.33 (dashed- 
dotted line), and x=0.40 (dotted line) in comparison with 
the corresponding experiments (inset) are plotted in Fig. 
3. It is shown that at low doping x=0.30 (solid line) and 
x=0.33 (dashed-dotted line), the spin-lattice relaxation 
rate 1/TiT reduces with decreasing temperature in the 
temperature region 0.2J sa 50K < T < 0.7J w 160K, 
which is self-consistent with the integrated dynamical 
spin susceptibility shown in Fig. 1 and Fig. 2, and is 
the result of the opening of normal state gap. However, 
in the low temperature region T < 0.2 J, the spin-lattice 



relaxation rate I/TiT follows a Curie- Weiss-like behav- 
ior due to the existence of strong AF fluctions. However, 
at the doping x=0.40 (dotted line), although the normal 
state gap still exits at high temperature. Curie- Weiss-like 
behavior at low temperature disappears, indicating that 
AF-fluctuations have become weak with the increasing 
doping. These results are in good agreement with the 



recent experiment; 



,11.12 



At the high doping x > 0.4, 



the neutron scattering studies and time-of-flight experi- 



ments on Nao.75Co02 and Nao.82CoO: 



9,10 



respectively. 



have observed the ferromagnetic correlation in the 2D 
C0O2 plane. Therefore, the magnetic fluctuation in the 
2D C0O2 plane has been affected by the electron doping, 
i.e., the magnetic fluctuation is doping dependent. The 
later ab-inito band structure calculations for Na2;Co02 
have given a certain critical value Xm = 0.56 ~ 0.68. 
The magnetic susceptibility in the C0O2 plane shows a 
tendency towards the AF state below Xm and ferromag- 
netic state above Xmr^. 



IV. SUMMARY 

In summary, we have studied the dynamical spin re- 
sponse in the electron doped cobaltates within the t-J 
model. We have shown that due to the presence of normal 
state gap, the particularly universal behaviors of inte- 
grated dynamical spin structure factor seen in the doped 
cuprates at the low energy, are absent in the doped Mott 
insulators on a triangular lattice. We have also shown 
that the spin-lattice relaxation rate I/Ti divided by T 
reduces with decreasing temperature in the temperature 
region above 0.2J«50K, then follows a Curie- Weiss-like 
behavior at the temperature less than 50K originating 
from the strong antiferromagnetic fluctuations, in quali- 
tative agreement with experiments. Our results also in- 
dicate that the electron doping plays a important role in 
altering the magnetic properties of the C0O2 plane. 
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